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I.  INTRODUCTION 


In  this  paper  a discussion  is  given  of  the  noise  sources  and 
sensitivity  limits  of  the  rf  biased  SQUID.  Both  the  low 
critical  current  "inductive"  and  the  higher  critical  current 
"hysteretic"  devices  are  considered.  The  intrinsic  SQUID 
noise  is  treated  by  considering  the  internal  flux  response  of 
a SQUID  ring  to  an  applied  rf  flux.  While  approximate,  this 
approach  yields  the  most  direct  picture  of  the  nature  of  the 
intrinsic  noise.  Circuit  optimization  procedures  for  maximiz- 
ing SQUID  performance  under  various  operating  conditions  are 
also  outlined.  Finally  the  effects  of  non-ideal  weak  link 
behavior  unon  SQUID  operation  are  treated. 


II.  SQUID  RING  RESPONSE 

A basic  understanding  of  the  origins  of  the  intrinsic  noise 
of  an  rf  SQUID  can  be  obtained  through  study  of  the  general 
nature  of  the  response  of  an  isolated  SQUID  ring  to  quasi- 
static and  time  varying  applied  magnetic  fields.  In  this 
section  wo  briefly  sketch  the  essential  details  of  this 
response,  emphasizing  those  features  which  ultimately  affect 
SQUID  sensitivity.  More  complete  discussions  can  be  found 
elsewhere  in  the  literature. 1-7  We  begin  by  considering  the 
case  of  very  low  frequency  or  dc  applied  flux.  The  general 
case  of  SQUID  ring  response  to  a combined  rf  and  dc  applied 
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flux  is  then  treated. 

II. 1.  Quasi-Static  Applied  Flux 

A standard  schematic  representation  of  an  rf  SQUID  ring  is 
shown  in  Fig.  1.  A superconducting  ring  of  inductance  Lg  is 
closed  by  a superconducting  Josephson  junction  or  weak  link. 
The  usual  shunted  junction  model  of  the  weak  link  is  assumed, 
where  a superconducting  element  with  a current-phase  relation 
i f (0)  is  shunted  by  a resistor  R and  capacitor  Cc.  For 
proper  SQUID  operation  it  is  necessary  that  the  weak  link  be 
overdamped  which  requires  Cg  < (4>0/2tt ic)/Rg2-  In  discus- 
sion that  follows  we  assume  this  is  the  case  and  that  hence 
C can  be  neglected  when  determining  SQUID  ring  response. 
Ideally  the  current-phase  relation  is  sinusoidal  although  in 
practice  this  is  often  not  the  case.^  Here  we  will  assume 
f(0)  = sin8;  the  effect  of  non-sinusoidal  f (0)  will  be  dis- 
cussed later. 

When  an  external  flux  $x  is  applied  to  the  ring,  the  response 
of  the  internal  flux  4 is  given  by 

1)  Lg$/Rg  + ♦ s $x  ~ Lsicsin(27f<J>/<}>0) 

Q 

As  discussed  by  Kurkijarvi  Eq.  (1)  describes  the  motion  of  <p 
in  the  potential 

2)  U(4>,$x)  = 2/2Lg  - ($o/2mccos(2n*/<fr0) 

In  the  absence  of  thermal  fluctuations  and  in  the  limit 
A /a  <<  R /L  the  total  current  i flowing  in  the  ring  is  the 

X o s s 

equilibrium  supercurrcnt  ig.  Under  these  conditions  Eq.  (1) 
reduces  to 


3)  4 - *x  - Lgicsin(2int>/<|.0) 


FIGURE  1.  A schematic  representation  of  an  rf  SQUID  ring. 
Idoaliy  the  current-phase  relation  f (—  2 n ^/s»0)  is  sinusoidal 

and  C8  sufficiently  small  that  the  weak  link  is  very  over- 
damped. 


FIGURE  2.  Quasi-static  internal  flux  versus  applied  flux 
response  curves  for  various  values  of  critical  currents. 
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As  illustrated  in  Fig.  2,  Eq.  (3)  in  general  describes  a 
snake-like  curve  in  the  $ vs  $x  plane  which  becomes  reentrant 
for  Lic  > ♦0/2tt  . For  Lic  < $Q/2v  the  system  has  only  one 
potential  minimum  and  $ is  a continuous  function  of  $x.  For 
Lic  1 $0/2ir  the  system  has  multiple  fluxoid  quantum  states. 
Segments  of  $ vs  $x  curves  such  as  in  Fig.  2 with  positive 
slope  correspond  to  regions  of  local  potential  minima,  while 
those  segments  with  negative  slope  correspond  to  local  poten- 
tial maxima.  This  is  shown  most  clearly  in  Fig.  3 where  the 

potential  U(<J>,$  ) is  plotted  for  the  case  Li  - 2.0$  . If, 

X c o 

for  example,  we  initially  set  $ = 0 and  $ = -1$  as  in  Fig.  3 

x o 

and  then  increase  $x,  in  the  absence  of  thermal  fluctuations 

$ will  follow  the  potential  minimum  along  the  curve  given  by 

Eq.  (3)  until  the  potential  barrier  vanishes  and  the  system 

makes  a rapid  transition  in  time  ~L_/RC  to  the  next  lower 

energy  state.  The  system  will  remain  in  this  new  fluxoid 

quantum  state  until  $x  is  varied  to  the  point  where  a new 

lower  energy  state  is  accessible.  If  an  oscillatory  flux  is 

applied  to  the  ring  with  amplitude  ><p  , where  $ = Li  + $ / A 

xc  xc  c o 

and  is  the  critical  applied  flux,  the  system  will  follow  a 
discontinuous,  hysteretic  path  in  the  $ vs  $x  plane  as  indi- 
cated in  Fig.  4. 

II. 2.  rf  Applied  Flux 

In  general  when  the  SQUID  ring  is  to  bo  operated  as  a flux 
detector,  a high  frequency  biasing  flux  $ sincut  and  a "dc" 
signal  flux  $x  ’ is  applied  to  the  SQUID  ring.  It  is  the 
non-linear  response  of  the  SQUID  ring  to  this  applied  flux 

4>  *x  m ♦*»Bin“t  ♦ ♦xd°) 

which  ultimately  yields  the  SQUID  signal. 

This  response  can  be  determined  by  solving  the  equation  of 
motion  given  by  Eq.  (3).  An  approximate  solution  can  be 
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FIGURE  4.  An  example  of  the  hysteretic  path  followed  by  <p  in 

the  piano  for  a slowly  vnryinq  periodic  applied  flux.  Tn 

this  case  the  critical  current  l.sic  = 2<f>0. 


obtained  in  a straightforward  manner  by  assuming  an  internal 
flux  response 


5)  4 = $msin (wt-0)  + $ ^dc^ 


While  this  approximation  neglects  higher  order  terms  in  u>  it 
retains  the  essential  features  of  SQUID  ring  behavior. 
Assuming  this  response,  three  coupled  equations  are  obtained: 

6a)  *m[Ctn0  + "W  = *xm/sine  " 2ctn0Lsic 

XJ1  (2ir*m/*o)  cos  (2Tr<1>  (<3C)  /*o) 

ctn0  - R /uiL  [1  + 2L  i 
6b)  ss  sc 

XJ1  (2lr*n/*o}  cos  (2ir*  (d<=)  /*o)  1 
Vdc)  = *<dc>  - L i 

6c)  X s c 

xJo(2Tr^m/$0)sin(2Tr4>'  V^0) 


These  equations  do  not,  in  general,  have  simple  analytical 
solutions.  It  is  however  instructive  to  examine  Eqs.  (6)  in 
various  limits  and  to  compare  the  results  to  complete  digital 
solutions  of  Eq.  (1). 

In  the  usual  case  of  SQUID  operation  the  bias  frequency 
u <<  R_/L_.  This  paper  will  only  consider  this  low  frequency 

S S Q 

situation.  The  case  whore  w ~ R /L  is  discussed  elsewhere. 

s s 

When  h>  <<  R_/L_  the  internal  flux  response  to  # sinwt  follows 

S S «\ul 

the  quasi-static  <p  vs  <f>x  response  curve  as  given  by  Eq.  (3) 
and  as  shown,  for  example,  in  Fig..  2.  The  effect  of  the  dc 
flux  <^dc^  is  simply  to  shift  the  equilibrium  point  about 
which  the  internal  flux  is  modulated  at  frequency  u>.  In  dis- 
cussing this  response,  throe  separate  critical  current  regimes 
can  be  considered: 
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FIGURE  5.  The  peak  amplitude  of  the  internal  flux  response  <)>,„ 

to  an  applied  rf  flux  <frx  * ♦xmsinut.  The  response  function  is 
plotted  for  $£dc)  * o,  S<J>0  and  VJ>0- 


jL 
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1.  Li  <<  A /2n 
a c To 

In  this  very  low  critical  current  limit  Eqs.  (6)  yield  the 
solutions 


0 = 0 

*<dc>  5 *(dc> 

T Yx 

= ♦*»  - 

The  internal  flux  response  is  in  phase  with  the  applied  flux 
and  the  amplitude  of  is  weakly  modulated  by  : the  mag- 

nitude of  this  modulation  depending  linearly  upon  L i and 

periodically  upon  <b  . This  solution  was  first  derived  by 
10  . , . . 

Hansma  in  his  original  analysis  of  the  low  critical  current 
rf  SQUID. 

To  illustrate  the  nature  of  the  SQUID  ring  response  in  this 
low  critical  current  regime  the  peak  amplitude  of  the  internal 
flux  response  to  = 4>xn,sinuit  as  obtained  from  digital  solu- 
tions of  Eq.  (1)  is  plotted  in  Fig.  5 versus  <p  for  the  three 
cases  if>x  = 0,  h<t>Q,  h <+>0 - This  result  was  obtained  assuming 
Lic  = 0.024>o  and  <*>Ls/Rs  = 0.1.  While  the  peak  amplitude  of 
the  internal  flux  has  both  a slightly  different  periodicity 
than  docs  $ , as  given  by  Eq.  (7c) , and  a scmewha-  large 
amplitude  at  high  <J>xm  drive  levels.  Fig.  5 does  indicate  the 
basic  nature  of  the  ring  response  for  very  small  critical 
current. 

2.  L i * ♦ /2ir 
s c o' 

As  Lgic  -*•  <PQ/2v  screening  currents  in  the  SQUID  ring  become 
important  and  Eqs.  (7)  no  longer  satisfactorily  explain  the 
SQUID  ring  response.  Recently  Soeronscn, 11  Erne  et  al.*2  and 
Hansma  et  al.A^  have  all  investigated  this  regime  through 
various  analytical  means  and  have  been  successful  in  obtaining 
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FIGURE  6.  the  peak  amplitude  of  the  rf  response  to  the  applied 
rf  flux  ,(,x  = <(.xn,sinwt  for  (a)  Lgic  = 0 - 1^0 , (b)  L..ic  = 0.2;o, 

(c)  Lsic  = and  (d)  L.jic  *-  0.5(|>o.  Due  to  the  relatively 

high  drive  frequency  u>I,s/Ra  l 0.1,  the  discontinuous  regime 


actually  begins  at  Lsic  s 0.2<|>Q  rather  than  at  Lsic  = <t>Q/2v . 
The  curves  are  plotted  for  - 0,  *s<{>0 , !s#0. 
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quite  accurate  descriptions  of  the  SQUID  ring  behavior.  For 
our  purposes  here,  however,  it  is  sufficient  simply  to  discuss 
the  general  aspects  of  the  $ vs  $x  response  as  deduced  from 
Eq.  (6)  or  from  digital  solutions  to  Eq.  (1). 

It  is  seen  from  Eq.  C6b)  that  for  sufficiently  small  u>L_/R_  <J> 
remains  essentially  in  phase  with  4>x  for  all  Lgic  < <t>0/2ir. 

The  amplitude  of  the  response  at  frequency  u can  be  ob- 
tained from  numerical  evaluations  of  Eqs.  (6)..  Alternatively, 
the  peak  amplitude  of  this  internal  flux  response  to 

t m ♦ sinwt  can  be  calculated  from  Eq.  (1) . Results  of  such 

x xm 

calculations  are  shown  for  several  values  of  Li  in  Fig.  6. 

s c 

Again  the  distinction  between  the  p aak  response  and  the 
component  at  frequency  to  is  slight:  a somewhat  different 
periodicity  and  a diminishing  amplitude  at  higher  drive  levels 

3*  Vc  > V2ir 

Once  Lgic  exceeds  4>0/2ir,  the  hysteretic  nature  of  the  quasi- 
static internal  flux  response  to  a periodic  applied  flux 
results  in  <J>m,  as  given  by  Eq.  (6}  , becoming  a discontinuous 
function  of  4*^-  (This  effect  occurs  at  progressively  larger 

values  of  L_i  as  <oL  /R_  1.  In  Fig.  6 where  toL  /R  = 0.1, 

s c s s s s ^ 

the  discontinuous  response  commences  at  Lsic/$0  : 0.2.  ) In 
this  discontinuous  regime  it  is  most  appropriate  not  to  pursue 
further  the  analysis  of  the  internal  flux  response  from  the 
point  of  view  of  a sinusoidal  oscillation  at  frequency  <d  but 
rather  to  rely  on  digital  solutions  of  the  peak  amplitude  of 
the  internal  flux  oscillation  as  a function  of  4^.  Examples 
of  this  peak  response  are  given  in  Fig.  6. 

II. 3.  Effect  of  dc  Applied  Flux 

A major  effect  of  the  dc  applied  flux  upon  the  <P~<PX  rf  re- 
sponse curve  is  the  periodic  modulation  of  the  values  of  $ 

for  which  d4  /d<p is  a maximum  (see  Fig.  6).  This  is  the 

ru  xm 
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case  regardless  of  the  value  of  Li  and  regardless  of  whether 
the  peak  amplitude  of  the  response  is  detected  on  the  compo- 
nent at  frequency  <u.  In  general,  in  an  rf  SQUID  system  the 
SQUID  ring  is  coupled  to  a biasing  circuit  in  such  a way  so 
that  those  regions  of  the  response  curve  where  d^/d^xm 

is  a maximum  present  maximum  opposition  to  changes  in  the  flux 
applied  by  the  biasing  circuit,  resulting  in  the  "steps"  in 
the  SQUID  circuit  rf  I-V  characteristic.  It  is  the  periodic 
modulation  by  <J>xdc^  of  the  location  max  of  these  regions 
of  maximum  opposition  that  yields  the  signal  in  the  rf  SQUID. 
This  is  discussed  in  somewhat  more  detail  in  section  IV. 


Thus  the  signal  that  can  be  obtained  from  the  SQUID  ring  de- 
pends directly  upon  the  amplitude  of  the  modulation  of  the 
response  curve  about  for  constant  <f>m-  In  the  very  low 

critical  current  regime,  L i <<  $ /2ir,  the  maximum  value  of 
this  modulation  amplitude  for  a small  s .gnal  flux  6 <J> ' ' can 
be  obtained  directly  from  Eq.  (7).  Assuming  that  the  ring  is 
biased  near  a point  where  (2ir$xn/$0)  *-s  a maximum  and  that 
^(dc)  kas  been  get  near  the  value,  = (n+%)$  , for  which 

JC  aC  O 

maximum  small  signal  sensitivity  is  obtained,  then  for  a small 
signal  6$xdc^  the  modulation  is 


8)  6 4 (6<t(dc))  = 

1 yxm,maxv  Tx  ' 


where  J,  is  the  first  maximum  of  the  first  order  Bessel 
i,m 

function  J^. 

Por  larger  critical  currents  the  maximum  modulation  amplitude 
about  ^ is  more  difficult  to  obtain  exactly.  However  a 
value  which  is  correct  to  within  a factor  of  order  unity  can 
be  obtained  by  employing  the  approximation  which  is  shown 
graphically  in  Fig.  7.  In  this  situation  the  rf  response 

curve  is  approximated  by  straight  line  segments  whose  slopes 
are  equal  to  the  local  minima  and  maxima  of  d^/d#  . it 

is  the  peak  value  of  the  rf  flux,  such  as  shown  in  Fig.  7, 


FIGURE  7.  The  straight  line  approximation  to  the  rf  re- 

sponse curve.  As  Lsic  ► $0/2n  the  accuracy  of  the  approxima- 
tion improves. 
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which  is  of  interest,  then  the  values  of  these  slopes  are 
(1  + 2irL  i /*  )‘*  and  (1  - 2nL  i _/4>_)  " 1 respectively.  If  the 

SCO  SCO 

response  component  at  frequency  u>  is  being  considered  then  the 
slopes  of  the  first  two  straight  line  segments  are 
[1  + 2LsicJ1'  (0) ]" * and  (1  - 2LsicJ1’ (3.7) ] * 1 . The  differ- 
ences between  the  two  cases  are  slight  and  in  the  discussion 
that  follows  only  the  peak  rf  response  function  will  be  con- 
sidered. The  modulation  of  the  response  curve  by  about 

6 is  also  assumed  to  be  triangular  which  is  only  strictly 
xro 

correct  for  Lsic  > $0/2ir. 


With  these  approximations  the  maximum  amplitude  of  the  modu- 
lation of  the  rf  response  curve  is  simply  wL_i_  for 

s c 

<<  1) . 
c)  < v2 

we  have 


2irL_i  5 and  <fr  /2  for  2ttL  i„  > <J>  (assuming  u>L_/R_ 

SCO  O SCO  o , 

Since  the  modulation  is  assumed  to  be  linear  for  <$<)>^ 


With  this  result  the  limiting  sensitivity  of  the  rf  SQUID  can 
be  determined  in  a straightforward  manner. 


III.  INTRINSIC  NOISE 

In  the  context  of  the  shunted  junction  model,  the  intrinsic 

noise  of  the  SQUID  originates  with  the  current  fluctuations  in 

the  weak  link  shunt  resistor  Rg.  When  the  SQUID  ring  is 

biased  at  frequency  u>  it  is  simply  the  ratio  of  its  response 

(dc) 

to  these  current  fluctuations  and  its  response  to  that 

determines  the  maximum  sensitivity  of  the  rf  SQUID.  Formally, 
the  effect  of  thc-so  fluctuations  can  be  determined  by  adding  a 

noise  current  term  L iM  to  Eq.  (1),  where  the  spectral  density 

8 M . |<  L 

of  this  noise  current  is  <i*(v)>  ■ • The  problem 

N U S S 
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then  becomes  one  of  calculating  the  fluctuations  of  $ in  the 
potential  U(4,4X)  for  4X  « ♦^sinwt  ♦ 4^dc* . 

XIX.  1.  Li  « 4/2* 
sc  o 

In  this  limit  and  for  wLg/Rg  <<  1 the  problem  reduces 

to  the  calculation  of  the  spectral  density  of  the  fluctuating 

flux  of  a classical  inductor  L . 

s 

10)  <«*a(v)%>  - <4kBT8L28/Rg)% 

As  mentioned  in  Section  II,  in  an  rf  SQUID  system  it  is  the 

modulation  of  the  location  of  the  rf  response  curve  about  the 

♦y,,  axis  by  that  provides  the  SQUID  signal.  As  a result 

it  is  necessary  to  cast  the  internal  flux  fluctuations  into 

the  form  of  equivalent  applied  flux  fluctuations  in  order  to 

obtain  the  intrinsic  noise  limited  flux  sensitivity.  This  can 

be  achieved  simply  by  employing  the  response  function. 

In  the  very  low  critical  current  regime  4_  z 4 for  all  <>  , 

m xm  m 

to  zeroth  order  in  2irLgi^4  • Thus  the  fluctuations  of  the 
SQUID  ring  are  equivalent  to  an  uncertainty  in  4^  of  spectral 
density. 

11)  <«4^ai(v)>,s  - <4kBT8/Ra)* 


Combining  Eq.  (10)  with  either  Eq.  (8)  or  Eq.  (9) , depending 
on  whether  the  SQUID  system  is  sensitive  to  the  SQUID  ring 
response  at  frequency  u or  to  the  peak  response,  the  maximum 
SQUID  sensitivity  or  the  minimum  observable  flux  signal  is  . 
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where  Af  is  the  post-detection  bandwidth  of  the  SQUID  system. 
This  result  was  originally  obtained  by  Danilov  and  Likharev14 
in  their  analysis  of  the  limiting  characteristics  of  the  low 
critical  current  SQUID. 

III. 2.  Lsic  S ♦0/2w 

For  larger  values  of  critical  current  the  fluctuation  ampli- 
tude of  ^ is  no  longer  independent  of  4 and  +x.  In  the  usual 
case  of  L„knT  <<  1 these  fluctuations  can  still  be  treated 
however  as  small  excursions  about  the  minimum  of  the  potential 
Therefore  the  spectral  density  of  the  fluctuations 
in  4 is  given  by 

13)  ■ Lg[l  + (2irLgic/$0) cos2ir$/<J>oJ ' * (4kBT8/Rg) ^ , 

a quantity  that  depends  directly  on  the  internal  flux  ♦ and 

hence  on  the  applied  flux  $x.  As  a result  the  uncertainty  in 

the  amplitude  of  the  internal  flux  response  depends  on  tho 

position  of  the  SQUID  ring  on  the  ♦m"$xin  response  curve.  This 

can  be  clearly  seen  in  Fig.  8 where  response  curves  as 

calculated  numerically  from  Eq.  (1),  with  a noise  current  term 

included,  are  shown  for  various  values  of  Li.  In  all  cases 

s c 

the  external  flux  amplitude  is  increased  sufficiently 

rapidly  in  the  calculation  so  that  fluctuations  in  are  not 
averaged  out  on  the  scale  of  the  drawing.  The  correlation 
between  the  fluctuation  amplitude  and  the  slope  of  the  re- 
sponse curve  is  apparent.  In  general  the  spectral  density  of 
the  fluctuations  in  $ is  given,  at  least  approximately  by 

14)  ‘l^(vl»l|SL,(d)1/d«XBll«tBVll1)11 

Since  during  rf  SQUID  operation  the  ring  is  biased  on  a region 
of  maximum  slope  on  the  ♦m“^xn  response  curve,  it  is  the  maxi- 
mum amplitude  of  this  spectral  density  that  stust  be  considered 
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I 


PICORE  8.  ♦ra-fXJIl  rf  response  curve  calculated  numerically  from 
a noise  term  included  vor  various  values  of  L.ic.  Again  the 
drive  frequency  u = 0.1RS/LS.  For  clarity  only  the  curve  for 
* n+0  is  shown  in  (a) . The  simulation  assumes  Ts  « IK 
and  L#  - 10" *H. 
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15)  >'*  = Ls(l  - ^'<Leic/*0)-'^*kBTs/Rs),, 


To  determine  the  maximum  SQUID  sensitivity  in  this  larger 
critical  current  regime  this  fluctuation  spectral  density 
again  mu3t  be  converted  into  equivalent  fluctuations  in  4> 

xm 

For  L_i  not  too  near  to  /2  and  for  L k_T  sufficiently 
SC  O S lJ  s 

small  so  that  noise  rounding  of  the  time-averaged  A -<f> 

m xm 

response  is  not  too  severe  this  can  be  accomplished  with  the 
use  of  the  noise  free  response  curve.  The  result  is 


16) 


= «i«VA|S 


Thus  the  equivalent  fluctuation  amplitude  in  $ remains  un- 


changed as  Li 
s o 


♦Q/2  , showing  that  the  simple  approach  of 


treating  the  noise  current  in  R as  the  source  of  an  applied 

noise  flux  6$^  is  valid  for  quite  large  critical  currents. 

Due  to  this  constancy  in  5<J>  (v)  the  minimum  detectable  flux 

signal  remains  as  given  by  Eq.  (12)  and  decreases  as 

Lgic  -*■  4>0/2tt.  While  this  result  is  no  longer  valid  once  Lgic 

is  sufficiently  close  to  <pQ/2ir , for  typical  values  of  SQUID 

parameters,  L_  - 10' ’H,  Td  - 4.2K  and  to  <<  L /R  , it  is  at 
s s s s 

least  approximately  correct  for  Lsic  S 4>o/4tt.  Consequently 

the  classical  noise  limit  6<j)^c^  s (4k_,T  L Af/R  ) **  can  be 

x B s s s 

rather  closely  approached  in  the  inductive  or  low  critical 
current  SQUID. 


III. 3. 


Vc  2 V2’ 


As  L8ic  approaches  from  below  the  value  at  which  the 

response  curve  becomes  discontinuous  (Li  = <p  /2v  for 

sc  o 

u»L  /R  <<  1)  the  amplitude  of  the  internal  flux  fluctuations, 
when  the  ring  is  biased  on  a region  of  maximum  slope  of  the 
response  curve,  eventually  increases  sufficiently  that  the 
intrinsic  noise  calculation  can  no  longer  be  treated  as  a 
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problem  of  classical  fluctuations.  Instead  it  must  be  treated 
as  a problem  of  random  transitions  between  quantum  states. 

This  effect  can  be  seen  in  Fig.  8 and  it  commences  somewhere 
below  Lsic  = $o/2n  at  a point  dependent  on  (LgkgT^)  . 

When  the  SQUID  ring  is  biased  at  or  near  the  region  of  maxi- 
mum slope  of  the  ♦m”,frxin  response  curve,  at  the  peak  of  the  rf 
cycle  the  internal  flux  <p  fluctuates  back  and  forth  between 
the  two  accessible  states  at  a frequency  which  is  much  greater 
than  the  bias  frequency  u),  assuming  ui  <<  R /L  . This  results 
in  an  uncertainty  in  the  peak  response  $ which  depends  on  the 
fluctuation  frequency.  As  Lgic  exceeds  <J>0/2tt  the  probability 
of  reverse  transition  diminishes  and  eventually  at  most  only 
one  forward  transition  in  $ can  be  made  at  the  peak  of  every 
half  rf  cycle.  As  a result  of  this  pronounced  decrease  in 
the  fluctuation  frequency  the  minimum  detectable  flux  of  the 
rf  SQUID  increases  until  the  maximum  SQUID  sensitivity  reaches 
the  value  calculated  by  Kurkijarvi  and  Webb^  for  Lgic  >>  <t>Q/2n. 


As  an  alternative  to  the  exact  calculation  of  Kurkijarvi  and 
Webb  the  minimum  detectable  flux  in  the  large  critical  current 
regime  can  be  discussed  along  lines  similar  to  those  used  here 
for  the  case  L i < tp/2ir.  Assuming  L i is  sufficiently 
large  and  <p ' 1 = n4>Q,  the  peak  amplitude  <J>m  of  the  rf  re- 
sponse to  <J>x  is  either  -0  or  <j>Q  during  any  rf  cycle,  depending 
on  whether  or  not  a transition  occurs  during  the  peak  of  that 
cycle.  If  the  probability  that  such  a transition  occurs  every 
cycle  is  p,  where  p depends  on  the  value  of  $ , then  the 

average  value  of  the  internal  flux  response  is  simply 


11)  «pn>  = p4,( 


and  for  frequencies  much  less  than  u>  the  spectral  density  of 

the  fluctuations  in  <t>  i3 

m 


18)  <6^ni(v)>}i  = (p(l-p)  |%0(2uA.)li 
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APPLIED  FLUX  AMPLITUDE  (<£xm/<£0) 


FIGURE  9.  A numerically  calculated  time  avcr.ujed  rf  response 
curve  lor  the  case  where  T = 0 and  T --  1.0K.  The  slope  ol  Liu 
dotted  curve  in  the  T ~ 1.0K  case  is  equal  to  the  fluxoid 
transition  uncertainty  as  calculated  by  Kurkijarvi  and  Webb. 


w 


To  convert  these  fluctuations  in  <J>m  into  equivalent  fluctua- 
tions in  $xm  it  is  necessary  to  use  the  response  curve 

which  has  been  time-avc raged  over  a period  >>  1/w.  An  example 
of  such  a time  averaged  response  curve  is  shown  in  Fig.  9. 

The  effect  of  the  averaging  is  to  remove  the  discontinuity 
from  the  response  curve,  replacing  it  with  a curve  of  finite 

slope  which  is  centered  at  a point  on  the  4>  axis  somewhat  to 

xm 

the  left  of  the  discontinuity  in  the  noise-free  case.  (In 

g 

general  the  amount  of  shift  depends  upon  i , T and  u).  ) This 

5 o 

curve,  which  is  basically  a straight  line  segment  over  most  of 
its  length  represents  the  probability  that  a fluxoid  transi- 
tion will  occur  once  every  cycle  for  a particular  value  of 
4>xm*  The  slope  of  the  curve  d<t)m/d<}>xro  is  then  simply  the  in- 
verse of  the  width  of  the  fluxoid  transition  uncertainty  as 

g 

calculated  by  Kurkijarvi.  Thus  the  equivalent  spectral 
density  of  fluctuations  in  where  the  SQUID  ring  is  biased 
on  the  region  of  maximum  slope  of  the  time-averaged  response 
curve  is 

19)  <64>2  (v)>*i  2 [p  ( 1-p)  ] **  (2ti/w)  **0 

xm 

where  a = Vc^bV^oV  ^ 

This  approximate  result  is  of  course  only  valid  for  biasing 
near  the  center  of  the  region  of  maximum  slope  of  the  response 
curve  whore  p ~ 0.9.  In  this  case  the  limiting  sensitivity  i ; 


20) 

which 


6$(dc*  i 0.5 
x,min 


'2’'Ws 

0)  J sc 

^o^c 

?/3 


is  essentially  the  result  of  K-W. 


9 


The  limiting  intrinsic  noise  for  an  rf  SQUID  is  thus  a minimum 

for  i slightly  less  than  the  value  required  for  the 

response  curve  l:o  become  di  scout  inuous,  which  is  Sn\.  i ./  0 

tor  wL  /R  <<  1.  As  i increases  above  this  value  the 
s o c 


1 
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intrinsic  noise  increases  rather  rapidly  toward  the  limit  cal- 
culated by  K-W  for  Li  >>  <j>  /2tt.  This  is  not  the  case  how- 
ever  when  the  bias  frequency  is  increased  so  that  toLg/Rs  > 1. 
For  wL  /R  - 1 it  is  straightforward  to  show  that  the  intrin- 
sic  noise  for  the  case  of  both  a continuous  and  discontinuous 
response  curve  can  very  closely  approach  the  limiting  classi- 
cal value.  Therefore  for  low  frequency  operation  uiLs/Rs  <<  1, 
if  the  intrinsic  noise  is  of  practical  importance,  which  to 
date  is  seldom  the  case,  then  ic  should  be  adjusted  so  that 
Lgic,  S <J>o/2n.  For  higher  frequencies  this  constraint  is 
largely  removed. 


IV.  CIRCUIT  CONSIDERATIONS 


In  the  usual  case,  the  rf  SQUID  circuit  consists  of  a current- 

driven  nigh  Q tank  circuit  of  inductance  and  capacitance 

Cfc.  The  circuit  is  inductively  coupled  to  the  SQUID  ring 

through  mutual  inductance  M and  the  circuit  is  driven  at  or 

near  its  resonant  frequency  ojq.  When  suitably  biased  and 

tuned,  the  tank  circuit  voltage  is  periodically  modulated  by 

6<p^ac^  , yielding  the  SQUID  signal.  A schematic  representation 

of  such  a SQUIT  circuit  is  shown  in  Fig.  10.  The  details  of 

the  tank  circuit  response  in  both  the  low  critical  current^ 

4 6 9 

ana  high  critical  current  ' ' regimes  have  been  extensively 
discussed  in  the  literature.  Hero  only  a few  basic  tea  tv.  s 
will  bo  sketched  so  as  to  establish  a qualitative  picture'  of 
the  various  noise  sources  and  optimum  coupling  conditions  for 
tne  rf  SQUID 


IV. 1.  rf  I-V  Characteristic 


In  general  terms,  the  rf  tank  circuit  is  driven  at  a frequency 
<»)  sucli  that  changes  of  the  rf  flux  amplitude  within  the  SQUID 
ring  present  maximum  opposition  to  change  in  the  applied 

1 lux.  As  d i seusned  by  ll.msm.i  and  ''  in  I ho  very  low 

critical  current  regime  where  the  flux  response  it  lrcqu.ney 
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is  in  phase  with  the  drive  flux  this  requires  that  the  bias- 
ing frequencies  to  be  cither  of  the  half  power  frequencies  of 
the  resonant  circuit.  For  larger  values  of  ic  when  screening 
currents  in  the  SQUID  ring  become  important  the  circuit  re- 
sponse is  no  longer  anti-symmetric  about  o>  and  maximum  signal 
is  obtained  by  tuning  to  the  larger  half  power  frequency.  As 

L i -*•  <p  /2ii  the  tuning  becomes  less  critical  and  for  L i 
s c o SC 

>>  <J>  / 2 ii  the  tank  circuit  which  is  now  sensitive  to  the  peak 
<f>m-$xm  resPonse  tunes  broadly  about  its  resonance  frequency  uiQ. 


It  is  a relatively  straightforward  procedure  to  construct  an 
approximate  rf  I-V  characteristic  of  an  optimally  tuned  SQUID 
system  by  using  the  simple  linear  approximation  to  the  time- 
averaged  4>m"4>xm  response  curve  such  as  shown  in  Fig.  7.  The 
basis  for  this  construction  is  the  following: 


When  the  tank  circuit  current  amplitude  Ifc  is  such  that  the 
rf  flux  applied  to  the  SQUID  ring,  <j>xm  = MIfc,  is  brought  to 
the  point  where  maximum  slope  is  found  in  the  d>  -d  response 
curve,  an  additional  impedance  is  brought  into  play  throu  h 
the  induction  by  the  SQUID  ring  of  a current  in  the  tank 
circuit.  This  current  is  directly  counter  tv")  the  current  in- 
duced by  the  external  drive  current  J,^.  As  a result  a gi  eater 
rf  drive  level  is  required  to  move  <&  and  hence  I.  through 
the  region  of  steepest  slope  in  the  <5  - response  curve  than 
through  the  region  of  minimum  slope.  Consequently  a "step"  is 
created  in  the  rf  I-V  characteristic  since  the  tank  circuit 
voltage  V,  = wL  = wL  $ /M.  The  length  of  this  step  in  tor  s 

u L L AM 

of  the  external  drive  current  is,  to  a fair  approximation, 


21) 


D , step 


L.L  , 
t s 


A.J> 


m 


a ib 

MQ  Fxm 


where  A<J>m  and  A^xm  are  the  extent  in  4>m  and  ^xrl*  respectively, 

covered  by  I he  steeper  straight  line  segments  ol  l he  i ■ *.»><•»?'  . 
iMivi1.  (S.e  Fig.  /.)  In  ( ! »<  ■ b w <vii  ical  cumin  i.  , i tu<  , o; 

in  Fig.  7,  we  have  A^/A^  = (1  - 2nLsic/<J>0)  ' 1 


, although  to  be 
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accurate  the  harmonic  response  function  should  be  employed 

for  which  A<J>m/A<J>xm  = [ (1  - 2LgicJ'1  ^)  (2Tr4>m/4>0)  ] * 1 gives  the 

slope  of  the  first  straight  line  segment  of  the  approximate 

response  curve.  For  L i > <J>  / 2n,  A<>  /A4>  is  the  reduced 

sc  o xm  xm 

fluxoid  transition  uncertainty  a/< p . The  rf  voltage  change 
along  the  step  is,  of  course, 


wL. 

22)  V.  . = A4> 

t,  step  M Txm 


The  various  steps  in  the  rf  I-V  characteristic  are  separated 
by  "risers".  The  length  of  these  risers  in  ID  can  be  esti- 
mated from  the  fact  that  from  the  end  of  one  full  step  to  the 
beginning  of  another  ID  must  increase  by  a minimum  amount 


23)  I - — (4>  - A<1>  ) 

D, riser  MQ  yo  Hxm 


while  the  tank  circuit  voltage  changes  by  an  amount 

w L*_ 


24)  V, 


t , riser 


™ (<f>  -A0  ) 

M o rxm 


Hy  determining  the  changes  in  position  of  the  various  steps 
and  risers  on  the  rf  I-V  characteristic  as  $xdcd  varies,  it 
is  quickly  found  that  the  voltage  modulation  by  the  applied 
flux  signal  of  the  positions  of  the  rf  steps  at  con- 
stant rf  drive  level  is 


<5V.  (6<J>  *ic))  = ~~  (1  - A0  /A -t>  ) <S$(dc)  , H(dC)  <•  0 /2 

t rx  M xnv  m x ' x o 


or 


25) 


6Vt  = M 


wLt|2,rLsic1 


603Jcic)  , 2,rLsic/4>0  < 0.5 


ml.. 


oV.  5 — — - (1  - U/v>  > I 


(dc) 


t ' M '*  ^ o'  ”rX  ' 2’V«/+o 


Win.  1 e Kqs.  (20)-(25)  are  clearly  only  approximate  , par  l ieu- 
Jarly  for  the  low  critical  current  case,  they  can  bo  used  to 
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construct  rf  I-V  characteristics  which  agree  reasonably  well 
with  I-V's  which  we  have  obtained  from  complete  digital  simu- 
lations of  rf  SQUIDs.  Examples  of  such  I-V's  are  shown  in 
Figs.  11  and  12  for  both  the  low  critical  current  and  high 
critical  current  cases.  In  all  the  I-V's  shown  the  drive  fre- 
quency is  at  or  close  to  its  optimum  value.  In  general  it  is 
found  from  the  simulations  that  the  values  of  the  maximum  dc 
flux  modulation  amplitude  and  of  the  step  and  riser  lengths 
arc,  at_  least  for  the  first  step,  well  within  a factor  of  two 
of  the  values  given  above,  once  corrections  have  been  made  for 
the  relatively  large  bias  frequencies  that  were  necessary  in 
the  simulations. 

A rather  important  result,  that  can  be  extracted  from  either 
the  arguments  given  above  or  directly  from  the  simulations,  is 
that  the  inductive  SQUID,  where  27rL  i /d>  <1,  can  yield  a 

SCO 

signal  that  is  within  a factor  of  two  of  the  maximum  that  is 
obtainable  from  the  higher  critical  current  device.  Thus  if 
intrinsic  noise  does  become  a practical  problem,  it  might  be 
appropriate  for  the  critical  current  to  be  reduced  from  its 
usual  value,  L„i  - <b  , to  Li  ~ 0.5<{>  /2ir,  where  the  intrin- 
sic  noise  is  considerably  less,  assuming  u>  <<  Ug/JJc. . 

IV. 2.  Circuit  Noise  and  Optimum  Coupling 

In  general  there  are  three  independent  noise  sources  in  an  r(. 
SQUID  — the  rf  amplifier,  the  tank  circuit,  and  the  SQUID  ring 
itself.  KurkijarviJ  has  analyzed  the  problem  of  tank  circuit 
noise  for  the  case  of  the  large  critical  current  device. 

Using  the  arguments  given  in  the  previous  section,  his  results 
can  be  readily  generalized  to  cover  the  lower  critical  current 
case . 

Tno  essential  point  of  Kurkijarvi's  analysis  is  that  the  amp- 
litude of  the  tank  circuit  voltage  fluctuations  is  delei/nined 
by  the  bias  position  on  the  rf  1 -V  characteristic;  that  is. 


Irf  (arb.  units) 


Irf  (arb.  units) 


FIGURE  12.  Simulated  rf  I-V  characteristic  for  I»sic  = 0.5^o. 
In  (a)  T ■ 0;  (b)  T = 4.2K.  K?Q  * 1 and  the  voltage  scale  is 
the  same  as  in  Fig.  11. 
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the  origin  of  tank  circuit  noise  is  the  current  fluctuations 
in  the  equivalent  tank  circuit  shunt  resistor  Rfc.  The  re- 
sponse of  the  circuit  to  these  fluctuations  is  determined  by 
the  local  tank  circuit  impedance  dVfc/dID.  Since  during  SQUID 
operation  the  system  is  biased  on  an  rf  step,  the  tank  circuit 
voltage  fluctuation  spectral  density  is 


dV 


^(v)^  - dl 


t,step 
D, step 


Vt 

Qu)L. 


or 


26)  <6V*  (v)>**  = — (4kRT  QioL.)1* 

1 Ad>  + K2QA<t>  13  C 11 

yxm  v Mn 

where  K*  = M2/LgLfc  and  is  the  coupling  constant  between  the 

SQUID  ring  and  L.  . Now  as  2irLi  /d>  increases  from  zero  to 

greater  than  one,  A$xm  decreases  from  4>Q/2  to  its  limiting 

value  o/d>  while  A<f>  increases  from  4>  /2  to  -A.  Thus  in  the 
o m o o 

limit  of  large  critical  current,  Eq.  (26)  reduces  to  the 

7 16 

Kurkijarvi  result.  ' Provided  the  critical  current  is  not 
too  small,  2mL  i /<p  i % and  that  K*Q  2 l,  Eq.  (26)  can  be 

3 C O 

approximated  by 

. (Ad)  /A<1>  ) , 

27)  <6V2  (V)>*  * — — (4k  T QwL  ) ** 

T K2Q  B t t 

where  A4>xn/Hm  = l-2irLsic/d>o  for  small  ic  and  Ad>xjn/Ad'm  = o/d>Q 
£or  Lsic  >:>  */2ir* 


With  this  result  the  total  equivalent  noise  power  spectral 
density  for  the  rf  SQUID  can  be  written  as 

K*L_ 


<Vjl(v)> 


+ 8 

"’Lt  (1  ' (1  ‘ 
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where  <6V^(v)>^  = spectral  density  of  the  voltage  noise  in  the 
amplifier  and  where  the  intrinsic  flux  noise  spectral  density 
<6«j>^(v}>1*  is 


< ^ 4> ( v)  >**  = 


2nL  i 

4knT  L 

s c 

B s s 

♦ 

R 

o 

s 

Li  < <p  /2n 
s c o 


29) 


a 


Vc 


>>  4>o/2tt 


Since  the  effective  tank  circuit  noise  temperature  Tfc  is 

usually  greater  than,  and  certainly  never  less  than,  T , Eqs. 

s 

(28)  and  (29)  show  that  for  the  inductive  SQUID  with  coupling 

K2Q  - 1 the  tank  circuit  noise  term  will  invariably  be  larger 

than  the  intrinsic  noise  term  provided  co  <<  R /L  . Only  for 

s s 

the  impractical  case  of  tight  coupling  does  the  tank  circuit 
noise  term  approach  the  intrinsic  noise  level,  which  however 
is  quite  small.  For  the  larger  critical  current  case  the 
intrinsic  noise  is  greater  and  the  ranking  between  the  two 
terms  depends  on  whether  the  ratio  L kuT./<}>  K Q $ v/2 . 

In  most  practical  SQUID  systems  the  amp’ifier  noise  term 
dominates.  To  minimize  this  noise  the  coupling  constant  K is 
reduced  to  some  minimum  value.  Usually  this  value  is  set  Ly 
the  requirement  that,  to  obtain  a continuous  modulation  of  the 
step  voltage  by  fi<|>  , it  is  necessary  that  I„  ^ 

For  Lgic  - <P0/2v  this  requires  that  K2Q  £ 1,  to  within  the 
accuracy  of  the  approximations  employed  here.  For  Li 

SC 

~ 0.5<J>o/2ir  the  extensive  rounding  in  the  actual  rf  I-V  charac- 
teristic which  is  not  accounted  for  in  the  approximations  used 
here  require  K2Q  to  be  somewhat  greater  than  unity  for  maximum 
modulation  of  Vfc  by  a flux  change  6<t>^dc*  = ♦Q/2. 

The  constraint  K2Q  - 1 is  not  rigorous  and  it  is  possible  I a 
increase  the  SQUID  signal  for  small  flux  charges  <5^dc*  *•  <l>0/2 
by  decoupling  beyond  K7Q  = 1.  This  is  most  effective  in  the 
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larger  critical  current  regime.  The  range  of  6<$>xdc*  for  which 
the  tank  circuit  voltage  is  responsive  decreases  as  ~1/K2. 
Still,  some  increase  in  the  small  signal  response  beyond  the 
value 


{vt  = "<°VLs>',<1  - 
is  possible. 


In  some  cases  decoupling  of  the  SQUID  ring  even  to  the  point 
K2Q  = 1 may  not  be  desirable  as  the  tank  circuit  noise  in- 
creases as  (K2Q)'X.  The  ratio  of  the  amplifier  noise  term  to 
tank  circuit  noise  term  is 


30) 


<«V*  (v)>* 
<6V2,(v)>Js 


K2Q<6V^(v) >* 

(QwLtkBTt),S 

A 6 
xm 

Since  the  maximum  value  of  A$m/A<f>xm  is  $10  for  typical  SQUID 
parameters,  for  sufficiently  high  tank  circuit  impedance  the 
tank  circuit  noise  may  become  a factor  as  improvements  in  rf 
amplifiers  continue  to  be  made.  In  this  case  tighter  coupling 
than  K2Q  = 1 may  be  required  for  maximum  sensitivity. 


V.  WEAK  LINK  EFFECTS  AND  ULTIMATE  SENSITIVITY 


Only  the  case  of  a sinusoidal  current-phase  relation  for  the 
weak  link  has  been  considered.  In  practice,  weak  links  that 
are  employed  in  SQUIDs  often  have  distinctly  non-sinusoidal 
current-phase  relations.  A typical  non-sinusoidal  current- 
phase  relation  is  shown  in  Fig.  13.  - The  result  of  such  de- 
partures from  ideal  behavior  is  not  drastic  but  it  is  detri- 
mental to  SQUID  sensitivity. 

The  major  effect  of  non-sinusoidal  current-phase  relations  is 
to  reduce  the  maximum  signal  amplitude  that  can  be  extracted 
from  a SQUID  system.  As  the  critical  phase  angle  0c  of  the 
current-phase  relation  increases  beyond  the  optimum  v/2  value 
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F1GUKH  13.  The  measured  non-si.nusoi.dal  current  phase  relation 
of  a l|i  square  tin  microbridcjc. 
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the  effect  is  to  decrease  the  amplitude  of  A4>m  on  the  steeper 


parts  of  the  rf  <P~<PX  response  curve. 


In  the  low  critical  cur- 
rent regime  the  result  is  a corresponding  increase  in  the 
minimum  coupling  constraint  on  K2Q.  For  0c  ~ tt  the  SQUID  can 
only  be  operated  in  the  discontinuous  or  hysteretic  mode  which 
has  the  higher  intrinsic  noise  amplitude.  Even  then  the  mini- 


mum coupling  constraint  Ir-tcp  ~ 


I . acts  to  decrease  the 
riser  _ 

In  addition  it  has  been  found / that 


in 

17 


available  SQUID  signal, 
the  fluxoid  transition  uncertainty  also  increases  with  criti- 
cal phase  angle  0 . Besides  increasing  the  intrinsic  noise 

c 

amplitude  directly,  this  imposes  a greater  lower  bound  on  the 
magnitude  of  This  results  in  turn  in  larger  tank 

circuit  noise  and,  if  the  transition  uncertainty  a ~ 0.3 ip 
a direct  decrease  in  maximum  signal  strength.  Experimentally 
it  has  been  found  that  even  for  extremely  non-sinusoidal 
current-phase  relations  the  transition  uncertainty  a seldom 
increases  beyond  a factor  of  5 of  the  value  for  a sinusoidal 
current-phase  relation.  It  appears  that  local  inhomogeneitios 
in  the  larger  weak  links  which  have  the  non-sinusoidal  rela- 
tions restrict  the  fluctuations  in  the  weak  link.  Consequent- 
ly, practically  any  weak  link  which  has  sufficiently  small 
critical  current  ic  ' <<>0/Ls  can  be  use^  in  an  rf  SQUID  in  the 
discontinuous  mode.  While  the  performance  is  degraded  if 

0 >>  it/2  a usable  signal  can  in  fact  usually  be  obtained, 

c 

provided  the  link  is  underdamped  and  w is  not  too  large. 


Finally  it  is  interesting,  although  perhaps  academic  at 
present,  to  determine  what  the  ultimate  intrinsic  SQUID  sensi- 
tivity might  be.  It  has  been  shown  that  in  general  when  dis- 
cussing SQUID  sensitivity  the  appropriate  sensitivity  figure 

is  not  the  total  equivalent  flux  noise  spectral  density 
* 

but  rather  one  which  takes  into  account  the  effici- 
ency by  which  flux  signals  can  be  coupled  into  the  SQUID  ring. 
In  this  case  the  SQUID  sensitivity  can  be  expressed  in  terms 
of  an  energy  resolution  6e  where  6E(v)  = <6$2  ( v)  >/K ' L„  where 
K'  is  the  coupling  constant  between  a signal  coil  and  the 
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SQUID  ring.  Assuming  that  K*  can  be  set  to  =1  for  the  most 

favorable  case  of  L i ~ d>  /2v , the  result  is 

sc  To 

31)  6E(v)  = 4kBTsLS/Rs 


For  an  ideal  Josephson  j unction  R = itA(0)/2ei  where  A ( 0 ) 

s c 

= 1.76kgTc  and  is  the  T = 0 energy  gap  of  the  superconductor 
and  where  the  temperature  of  the  weak  link  is  assumed  to 
be  <<T  of  the  superconducting  material.  In  this  case  we 
have 


32) 


6E(v) 


4ltBT8Ls(2elc> 


but  L i 
s c 


V211 


_h_ 
4ire  ' 


33)  6E(v) 


4T_ 


1. 76nT. 


_h_ 

2tt 


so 


Thus  for  any  rf  SQUID  system  which  has  an  ideal  Josephron 
junction  for  the  weak  link,  fabricated  from  any  superconductor/ 
and  which  is  operated  at  T $ 0.5T  the  minimum  intrinsic 

C 

energy  resolution  is  simply  due  to  the  zero  point  oscillations 
as  set  by  the  uncertainty  principle: 

<6<{>*(v)> 

34)  ~ h ~ 6.6>«10'3‘‘  Joules/Hz 

s 


In  terms  of  magnetic  flux  for  a typical  ring  inductance 
Lg  ~ 10'  *11  this  corresponds  to  an  ultimate  flux  sensitivity 

35)  <6<^(v)  >**  Z 4x10'  7$0/ Alz  , 

which  is  somewhat  more  than  a factor  of  ten  lower  than  the 
best  sensitivity  that  has  ever  been  achieved  in  practice, 
to  the  knowledge  of  this  author.  It  would  be  rather  ir.prc.j- 
sive  if  this  limiting  noise  figure  is  in  fact  ever  closely 
approached. 
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